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Systematic analysis of deformations in neutron-rich A'' = 14 isotones was done based on the 
method of antisymmetrized molecular dynamics. The property of the shape coexistence in ^*Si, 
which is known to have the oblate ground state and the prolate excited states, was successfully 
described. The results suggest that the shape coexistence may occur also in neutron-rich A'' = 14 
nuclei as well as ^*Si. It was found that the oblate neutron shapes are favored because of the spin- 
orbit force in most of A*" = 14 isotones. Q moments and E2 transition strengths in the neutron-rich 
nuclei were discussed in relation to the intrinsic deformations, and a possible difference between the 
' proton and neutron deformations in ^*Ne was proposed. 
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CN ■ I. INTRODUCTION 

^: 

^ \ The shape coexistence in ^^Si has been studied for a long time. At the early stage, the oblate and prolate solutions 
• of Hartree-Fock calculations are associated with the ground 0"*" state and the Oj^ state at 6.691 MeV, respectively 
' [1]. Experimentally, the ground band is known to be oblate from the sign of the quadrupole moment of the state, 
04 , while the the band- head Oj^ in the prolate band has been identified by gamma transition measurements [2] . 

The coexistence of oblate and prolate states in ^®Si can be described by the shell gap at = 14 and Z = 14 of 
Nilsson orbits in the deformed system [3]. A variety of shapes is considered to arise from the nature of the shell 
structure in iV = 14 and/or Z = 14 systems, where a half of the sc?-shell is occupied. It is natural to expect that 
\^ \ shape coexistence phenomena may appear also in neutron-rich nuclei with A^ — 14. From the coexistence of oblate 
' and prolate states in ^*Si, one may expect that this trend of the oblate and prolate neutron structure might survive 
0^ in other A^ = 14 isotones. The primary interest in the present paper is on the deformations of the proton and neutron 
. densities in neutron excess A^ = 14 nuclei. How does the oblate neutron structure in the ground band of ^^Si change 
' with the decrease of proton number toward the neutron-rich region ? Whether or not the coexistence of oblate and 
prolate deformations occurs ? The neutron deformations should be sensitive to the proton structure, therefore, the 
shape coexistence of the neutron structure may depend on the proton number. If a oblate neutron shape is favored 
also in neutron-rich nucleus as well as in ^*Si, one may meet another question: is the proton shape consistent with 
the neutron one ? In this sense, ^^Ne is an attractive nucleus, where the proton shape would be different from the 
^ ' neutron shape, because the prolate shape of proton density is favored in Z = 10 systems as seen in ^"Ne. 
f5 In the study of proton and neutron deformations, the experimental data of electric moments such as B(E2) are 

useful to extract information about the intrinsic deformations. Recently, B{E2; 2+ 0+) in ^^Si has been measured 
by Coulomb excitation [4]. The measured B{E2) in ^''Si is almost as large as the experimental value of B{E2) in ^^Ne. 
The ratio B{E2-^^'&i) / B{E2;'^^^e)< 1 is much smaller than a naive expectation based on a collective model picture 
H \ that B{E2) is proportional to {N/Zf « 2. If these nuclei are mirror symmetric, the ratio B{E2-^^^i) / B{E2\'^^'S>ie)< 1 
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lead to a possible difference between proton and neutron deformations in ^"^Ne. 

The coexistence of the oblate and prolate solutions in ^*Si has been confirmed by many theoretical calculations such 
as Nilsson-Strutinsky calculations [5], Hartree-Fock-Bogoliubov [6], alpha-cluster model approaches [7,8]. Although 
mean field approaches are useful for systematic study of deformations, their applicability to very light nuclei is not 
obvious. In the light nuclear region, we would like to remind the readers of the importance of cluster aspect, which 
closely relates with the deformations. The cluster aspect has been suggested also in light neutron-rich nuclei as well 
as stable nuclei. For example, it is well known that ^°Ne has a ^^O-t-a cluster structure, while the development 
of clustering was suggested in neutron- rich B isotopes. The cluster aspect and the shape coexistence of even-even 
N = Z nuclei has been studied by a-cluster models [8,7]. 7a-cluster model calculations were applied to ^*Si, and 
succeeded to describe an exotic shape with Ds/i symmetry [7,9], which is associated with the K"^ = 5^ band observed 
in gamma transitions [2] . Although the properties of the oblate and prolate states were described by the 7Q;-cluster 
models [8,7], however, many 7a-cluster calculations failed to reproduce the order of the oblate and prolate solutions 
in ^^Si. Namely, in most of the a-cluster calculations, the prolate solutions are lower than the oblate solutions except 
for few calculations with Brink-Beoker forces [7]. For the oblate property of the ground state in ^*Si, it is important 
to incorporate the effect of spin-orbit force, which is omitted within the a-cluster models. Furthermore, the cluster 
models are not suitable for the systematic study of unstable nuclei, because they replies on the assumption of existence 
of cluster cores. 
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For the systematic structure study of light nuclei, one of the powerful approaches is the method of antisymmetrized 
molecular dynamics (AMD). The applicability of this method for unstable nuclei has been proved in many works 
[10 13]. In addition to description of deformations and cluster aspect in light nuclei, this method has an advantage 
that electric moments can be directly related with structure change based on the microscopic treatment of spin-parity 
projection without introducing effective charges. 

In this paper, we study the deformations of = 14 isotones while focusing on the coexistence of the oblate and 
prolate neutron structure. We pay attention to reproduction of the order of oblate and prolate solutions, which 
coexist in ^^Si, and discuss the effect of the spin-orbit force. We analyze the systematics of deformations in the 
neutron excess nuclei in relation to the obscrvablcs such as Q-moments and E2 transitions strength. This paper is 
organized as follows. In the next section, the formulation of AMD is briefly explained. We show the theoretical results 
and give comparisons with the experimental data in Sec. IV. In Sec.V, we analyze the instrinsic structure, and discuss 
the lowering mechanism of the oblate state within the AMD framework. Finally, a summary is given in Sec.VI. 



II. FORMULATION 



Here we briefly explain the formulations of the present calculations. Details of the formulation of AMD methods 
for nuclear structure studies are explained in Refs. [10,12]. The present calculations are basically same as those in 
Ref. [10]. 

The wave function of a system with a mass number A is written by a superposition of AMD wave functions $amd- 
An AMD wave function is given by a single Slator determinant of Gaussian wave packets as, 

$amd(Z) = -^A{ipi, IP2, </9a}, (1) 

where the i-th single-particle wave function is written as, 

V>i = -^XiXiTi, (2) 

(^x^(rj) a exp{-u{rj - -^f}, (3) 

Xi = {l+^i)X1 + il-i^)Xl■ (4) 

The iso-spin function is fixed to be up(proton) or down(ncutron), and the orientation of intrinsic spin is fixed to 
be 1/2 or —1/2 in the present calculations as done in Ref. [10]. The spatial part is represented by complex variational 
parameters, Xii, X2i, Xsi, which indicate the centers of Gaussian wave packets. 

In the AMD model, all the centers of single-nucleon Gaussians are treated independently as the complex variational 
parameters. Thus, this method is based completely on single nucleons and therefore it does not rely on the assumption 
of the existence of chister cores. In the sense that a single AMD wave function is written by a Slator determinant of 
Gaussians, the AMD method is regarded as an extended model of Bloch-Brink cluster model [14]. Here we note that 
the Gaussian center is expressed by the 'complex' parameter Xj which contains a real part and imaginary part. It 
means an extension of the model space, because the degrees of freedom are twice of the case that Gaussian centers 
are given by real values as in usual cluster models. If we ignore the effect of antisymmetrization, the position Dj 
and the momentum Kj of the i-th single-nucleon wave packets are expressed by the real and imaginary parts of Xj, 
respectively as, 
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ixjr|,/.x,) Re(Xi) 
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K. . = 2nV^lmiX,). (5) 

(?'xj<pxj 

In the nuclear structure study with the AMD, the imaginary parts of Z are essential to describe the rotation motion 
of the system. They are important to incorporate the effect of spin-orbit force and to describe high-spin states. 

We perform energy variation for a parity-eigen state, P^^amd = ^'amd' Projected from an AMD wave function 
by using the frictional cooling method [10]. We consider the AMD wave function obtained by the energy variation 
as the intrinsic state, and total-angular- momentum projection(P^^^) is performed after the variation to calculate the 
expectation values of operators such as energies and moments. In the present calculations, the parity projection is 
done before the variation, but the total-angular-momentum projection is performed after the variation. In many of 
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A'^ = 14 isotoncs, two local minimum solutions are found in the energy variation. In such cases, we diagonalize the 
Hamiltonian and norm matrices, {Pmk'^amj:>\^\Pmk"^j^u) ^^"^ iPMK'^AMD\PMK"^AiD) ^^^^ respect to the 
obtained intrinsic wave functions (*!?amd>^amd) ^^'^ the ii'-quantum {K',K"). After the diagonalization the ground 
and excited bands are obtained. 

III. INTERACTIONS 

The effective nuclear interactions adopted in the present work consist of the central force, the spin-orbit force and 

Coulomb force. We adopt MVl force [15] as the central force. The MVl force contains a zero-range three-body 
force in addition to the two-body interaction. The Bertlett and Heisenberg terms are chosen to be 6 = /i = 0. We 
use a parameter set, case 1 of MVl force with the Majorana parameter m = 0.62. Concerning the spin-orbit force, 
the same form of the two-range Gaussian as the GSRS force [16] is adopted. The strengths of the spin-orbit force, 
ui = —uii = uis = 900, and uis — 2800 MeV are used. 

IV. RESULTS 

We calculate the natural parity states of TV = 14 nuclei,i'^B, ^Oc, ^^Q, '^'^Ne and ^^Si, with the AMD method. The 
width parameters z/=0.145, 0.15, 0.150, 0.155, 0.15 are chosen for ^^B, ^Oc, 229, 24]\fg ^nd ^^Si, respectively, so as to 
minimize the energy of each nucleus. 

In these A*" = 14 nuclei, we find two local minimum solutions with oblate and prolate deformations of neutron 
density, except for the results of ^^C and ^^O with uis = 2800 MeV. It signifies the trend of the shape coexistence of 
neutron structure in the A'^ = 14 isotones. The feature of the shape coexistence originates in the nature oi N = 14 
neutron structure. 

After performing the total-angular-momentum projection and the diagonalization for the obtained oblate and 

prolate states, we obtain rotational bands, K'^ = 3/2 j~ and 3/2^ bands in -'^^B, and K"^ = 0^ and 0^ bands in ^^Ne 
and 2® Si from the two local minimum solutions. In the calculations of 2''C and 22 with the stronger spin-orbit force 
Uis = 2800 MeV, the shape coexistence phenomena disappear and only one intrinsic state is obtained by the energy 
variation in each nucleus. In the calculations with uis = 2800 MeV, the lowest states of 220 and '^'^'Ne have small 
prolate deformations of neutron density. As mentioned later, these states with 'small prolate deformations' should be 
classificid as the 'oblate' states in tlie analysis of the shape coexistence in iV = 14 systems, because their properties 
are similar to those of the oblate states of other = 14 nuclei. Therefore, we call these lowest states of 22 O and 2'*Ne 
as 'oblate' states in this paper. 

The energies of the band- head states are shown in Fig.l. In the results with uu = 900 MeV, we find that oblate 
and prolate bands coexist in those A^ = 14 nuclei. The oblate and prolate bands are almost degenerate in 20 C and 
220, while in ^^B, ^^Ne and ^^Si the prolate bands are lower than the oblate ones. With the stronger spin-orbit force 
uis = 2800 MeV, the 'oblate' solutions are the lowest in all of these nuclei. Namely, the ground bands have the 'oblate' 
neutron shapes . The prolate bands appear as excited bands in ^^B, '^'^Ne and 2®Si. We can not find prolate solutions 
in 20C and 22 in the calculation with = 2800 MeV. 

The ground state of 2^81 is experimentally known to be oblate. The order of the oblate and prolate bands is well 
reproduced by the calculations with uis = 2800 MeV. On the other hand, the calculations with uu = 900 MeV tend to 
relatively overestimate the energy of the oblate bands, and they fail to reproduce the order of the oblate and prolate 
solutions in 2^81. Comparing the results with uu = 900 MeV and uu = 2800 MeV, it is found that the relative energy 
of the oblate and prolate states in A/' = 14 system is very sensitive to the strength of the spin-orbit force. In order 
to describe the oblate feature of the ground state of "^^Si, the effect of the spin-orbit force is significant as follows. 
In comparison of the results for the stronger spin-orbit force(u;s = 2800 MeV) with those for weaker one(u;s = 900 
MeV) in Fig.l, we find that the energy gains in the oblate states are as large as about 10 MeV, while the energies 
of prolate states are almost unchanged with the increase of the strength of the spin-orbit force. As shown in table I, 
the absolute value of the spin-orbit term is much larger in the oblate states than in the prolate states. It indicates 
that the energy gain of the spin-orbit force in the oblate state is essential to describe the level structure of 2*Si. The 
stronger spin-orbit force Uis = 2800 is appropriate in the present framework for describing the order of the oblate and 
prolate levels in "^^Si. We discuss the lowering mechanism of the oblate states in section V. 

Next, we show the level scheme obtained witli uis = 2800 MeV in Fig. 2. 

In ^^B, we predict the coexistence of oblate and prolate bands. In the previous AMD calculations of B isotopes 
[10], the prolate shape in the ground state was predicted in ^^B. This is because the adopted spin-orbit force was too 
weak in the previous calculations. On the other hand, when we use the stronger spin-orbit force, which reproduces 
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FIG. 1. Calculated binding energies for the oblate(o) and prolate(p) states. The energies are of the band-head states, 



in '^B and the 0+ states in ^"C, 



'Si. The results with uu = 900 McV{uis = 2800 MeV) are shown 



J± = 3/2- 

in left(right). Although the 'oblate'('o') states in ^^O and ^*Ne have small prolate deformations of neutron density, these states 
are attributed to the oblate neutron structure in iV = 14 systems. The details are described in the text. 
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the order of the oblate and prolate bands in ^^Si, we obtain the ground state of ^^B with the oblate neutron shape 
while the prolately deformed state appears above the ground state. As shown later, the excited prolate band has 
^^Li+^He-like cluster structure. Although the 2-neutron separation energy of ^^B is very small as about 0.1 MeV, we 
expect that the prolate states may exist as resonances because such cluster states might be stable against neutron 

decays. 

In ^^O, we obtain the ground state but can not find other excited states in the present calculations. This is because 
of the double shell-closed feature of ^^O. Namely, the intrinsic state has mostly spherical shape among the N = 14 
nuclei due to the effects of the neutron sub-shell closure and the proton magic number Z = 8. 

In the results of ^^Ne, the lowest 'oblate' band in the theoretical rc;sults corresponds to the ground band of exper- 
imental data as discussed in the next section. On the other hand, the excited prolate bands are not experimentally 
identified yet. In the present calculations, a side-band = 2+ of the oblate ground is suggested in addition to the 
excited prolate band. 

In ^^Si, the shape coexistence of the oblate ground band and the prolate excited band has been known. The prolate 
band, which starts from the band-head 0^ state at 6.691 MeV, has been identified in gamma-transition measurements 
[2]. The O2 state at 4.979 MeV in the experimental data is considered to be a vibrational excitation in the oblate 
state. The level structure of the oblate and prolate bands are reproduced by the present calculations. The level 
spacing between oblate and prolate bands can be reproduced in good agreement by using a slightly stronger spin-orbit 
force uis = 3200 MeV which gives 7.3 MeV excitation energy for the band-head state of the prolate band. 

The theoretical results of the root-mean-square matter radii, E2 transition strengths, and electric and magnetic 
moments are shown in Fig. 3 and tables II and III in comparison with the experimental data. The present calculations 
reasonably agree to the experimental data. 
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TABLE I. Expectation values(MeV) of the spin-orbit force in the oblate and prolate states obtained with uis = 2800 MeV. 
Each value is for the band-head state obtained by spin-parity projection from an single AMD wave function which corresponds 
to 'oblate' or prolate solution. The superposition of the oblate and prolate solutions is not done. 
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FIG. 2. The level structure of "B, ^"C, ^^O, ^*Ne and ^*Si calculated with uu='M)^ MeV. 



5 



3.4 




td 2.6 

B 

2.4 

24Ne 28si 

FIG. 3. Root-mean-square matter radii. Open circles indicate the radii of the lowest states calculated with uu = 2800 MeV. 
The experimental data are shown by square points. The data for ^®B, ^"C, ^^O and ^*Ne are those derived from the interaction 
cross sections [17], and the radius of ^^Si is derived from the charge radius taken from Ref. [18] measured by electron scattering. 



TABLE II. E2 transition strengths. The results are calculated with u;s=2800 MeV. The theoretical B{E2) in ^*Si arc 
evaluated by assuming mirror symmetry with ^''Ne. The experimental data are taken from Refs. [2,4,19]. " Since the collective 
4+ 2+ transition in the = Oj band of ^®Si is distributed over two transitions into 2+(7.381) and 2+(7.417), the 
experimental value in this table is a sum of the strengths of those two transitions in Ref. [2]. 
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TABLE III. Electric quadrupole moments and magnetic dipole moments. 
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V. DISCUSSION 



In this section, we analyze the deformations of proton and neutron densities in the intrinsic states and discuss 
their effects on the observables such as E2 transitions. The following discussions are based on the calculations with 
uis=2800 MeV. 

A. Intrinsic deformations 

As mentioned before, two local minimum states arc obtained in each system except for and ^^O. The deformation 
parameters /?, 7 for neutron and proton densities in the intrinsic states are shown in Fig. 4. The definition of [5, 7 are 
given in Ref. [20]. Figure 5 shows the distribution of matter, proton and neutron densities. As seen in Figs. 4 and 
5, the lowest states of ^^B, and ^^Si have oblate neutron shapes. Although the lowest states of ^^O and ^"'Ne 
have small prolate deformations of neutron densities, we notice that the characteristics of the neutron structure in 
these states are rather similar to those in the oblate states of other A'' = 14 nuclei. One of the remarkable features 
of the 'oblate' neutron structure is the larger energy gain of the spin-orbit force than the prolate neutron structure. 
In fact, the energy gains of the spin-orbit force in these states of ^^O and ^^Ne are as large as those in the oblate 
states in -'^^B, ^"C and ^*Si as seen in Fig.l and table I. It is considered that the oblatcly deformed neutron structure 
is not rigid but is somehow soft to vary into the small prolate deformation in such systems as ^^O and ^''Ne, which 
have the spherical and prolate proton shapes, respectively. In other words, the oblate neutron deformations can be 
modified because of the inconsistency with the proton deformations. It should be stressed that the features of the 
shape coexistence of neutron structure in A^ = 14 nuclei varies depending on the proton number. 

In ^^B, the lowest state has the oblate neutron shape and a triaxial proton shape, while in the excited state the large 
prolate deformation with a cluster structure is enhanced. As seen in Fig.5, the developed *He-|-^^ Li-like clustering is 
suggested in the prolate excited band. 

In the lowest state of ^°C, both the proton and neutron shapes arc oblate. The reason for the absence of the prolate 
solution in is understood as follows. As known in ^^C, the proton shape is oblate in Z = 6 nuclei. Because of the 
inconsistency with the oblate proton deformation, the prolate neutron structure is energetically unfavored comparing 
with the oblate neutron structure. 

In ^^O, the neutron deformation is the smallest among these A^ = 14 isotones due to the spherical proton shape 
which originates in the p-shell closure effect. Another local minimum solution with the large prolate deformation of 
neutron structure does not appear as well as in ^"C. 

In ^^Ne, two local minimum solutions with different neutron structures are obtained. The neutron structure of the 
lowest state has a smaller prolate deformation than that of the excited state. This smaller neutron deformation in the 
lowest state can be attributed to the oblate neutron structure in A'' = 14 systems, because as mentioned above the 
sensitivity of the energy of this state to the strength of the spin-orbit force is quite similar to that of the oblate states in 
other N = 14 isotones. On the other hand, in the excited state, the neutron structure has a large prolate deformation 
and is similar to the prolate excited state of ^*Si. Therefore, the coexistence of two local minimum solutions in ^"'Ne 
is associated with the shape coexistence of the oblate and prolate states in ^^Si. One of the unique features in ^^Ne is 
that the proton deformations are prolate in both the groimd and excited states. It is because of the nature of Z = 10 
systems. Due to the effect of prolate proton deformation, the oblate neutron structure varies into the small prolate 
deformation. As a result, the deformation parameters /3p and /3„ for the proton and neutron densities are inconsistent 
with each other as /3p > /3„. This leads to a difference between proton and neutron deformations in the ground state 
of ^^Nc. The details of the different deformations are discussed later in relation to the E2 transition strengths. 

In the results of ^*Si, the oblate and prolate states coexist. As already mentioned, the oblate feature of the ground 
state is reproduced with the stronger spin-orbit iorce{uis = 2800 MeV). The shape coexistence of ^*Si has been studied 
by 7a-cluster models [7-9], where exotic shapes have been suggested in addition to the normal oblate and prolate 
deformations. For example, a pentagon shape composed of 7 a clusters with D^h symmetry [9] is attributed to the 
K'^ = 5~ rotational band which has been observed in gamma transition measurements [2]. In the present calculations, 
all the centers of 28 single-nucleon Gaussians are treated as independent variational parameters without assuming 
existence of any cluster cores. By analysing the positions Di(i = 1 ^ 28) of Gaussian centers, we find that the oblate 
and prolate states consist of 7 a-like clusters in the present results. Schematic figures for spatial configurations of 
the oblate and prolate states are shown in Fig. 6. Due to the 7a clustering, the pentagon shape appear in the oblate 
solution as seen in Fig.5. We give a comment on 'a' clusters in the present results. Since the imaginary parts of the 
single-nucleon Gaussian centers are non-zero in the oblate solution, the a-like clusters are dissociated a clusters which 
differ from the ideal a clusters written by (Os)^ harmonic-oscillator wave functions. The details are discussed later in 
the last part of this section. 



7 



neutron 




proton 




FIG. 4. Deformation parameters of the intrinsic states $amd(Z) obtained with uis = 2800 MeV. The deformation parameters 
f3n,'yn{f3p,jp) for the neutron(proton) densities are plotted in the upper(lower) panel. Triangles denote the deformations in the 
ground bands of ^^B, ^°C, ^^O, ^*Ne and ^®Si, while circles indicate those in the excited bands of ^^B*, ^''Ne* and ^*Si*. 
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FIG. 5. Density distribution of the intrinsic state #amd(Z). The axis of the intrinsic frame are chosen as {xx) < (yy) < (zz) 
and {xy) = (yz) — (zx) = 0. The intrinsic system is projected onto the zy-plane. The density is integrated along the transverse 
axis X. The densities for matter, protons and neutrons are displayed in the left, middle and right figures, respectively, (xx), 
{yy}, {zz) for matter, proton and neutron densities are written below the figures. The unit of the box frame size is 10 fmxlO 
fm. 



(a) (b) 




FIG. 6. Schematic figures for the positions(Di) of the Gaussian centers in the oblato(a) and prolate(b) states of ^*Si. The 
positions Di(i = 1 ~ 28) are distributed in 7 groups which associate with the 7a-cluster configurations. 
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B. i»B 



The structure of B isotopes has been studied with the AMD method in Ref. [10] . In the previous results, it was 

predicted that chister structure develops as a neutron number increases from "'^^B toward the neutron-drip line ^^B 
due to a large prolate deformation of neutron structure. In the present results, however, the ground state of ^^B 
has the oblate neutron structure, which is inconsistent with the previous results. In the previous calculations of ^^B, 
the oblate solution was not lower than the prolate solution because the adopted spin-orbit force uis=900 MeV was 
too weak. As explained before, the relative energy of the prolate and oblate solutions is sensitive to the strength of 
spin-orbit force. Therefore, it is important to use a proper strength of the spin-orbit force to describe the feature of 
shape coexistence in N = 14 systems. By using the stronger spin-orbit force U(s=2800 MeV, which can reproduces 
the level scheme of the oblate and prolate bands in ^^Si, the oblate band becomes lowest in ^^B while the prolate 
band appear as the excited band (Fig. 2). As a result, the neutron structure of the groimd state of ^'^B is predicted 
to be oblate in the present work. The oblate property of the ground band of ^^B is contrast to the features of other 
neutron-rich B isotopes(^^B and ^'^B), where prolate neutron structure is favored [10]. 

Because of the oblate neutron structure in ^''B, the proton deformation is smaller, and the clustering is weaker in 
the ground state than the excited prolate state as seen in Fig. 5. The deformation of the neutron structure is reflected 
in the electric quadrupole moment(Q) through the proton deformation. As shown in table III, the calculated Q{3/2~) 
for the oblate ground state is much smaller than that for the prolate excited state. The experimental measurement 
of the Q moment of ^^B is required to confirm the oblate neutron deformation. Contrary to the Q moments, the 
calculated magnetic dipole moments are not so much different between oblate and prolate states. 



Ne is an interesting nucleus, because the proton structure in Z = 10 system may favor a prolate shape as known 
in ^°Ne. Considering the coexistence of oblate and prolate shapes in A'' = 14 neutron structure, it is natural to 
expect that the neutron shape and the proton shape may compete. The present results suggest that the different 
kinds of neutron structure coexist in ^^Ne. One is the smaller prolate deformation in the ground band, and the other 
is the larger prolate deformation in the excited band. The former one is attributed to the oblate structure in the 
AT = 14 systems. We consider that the change of the neutron structure from the oblate shape into the small prolate 
deformation is understood as the oblate neutron structure is somehow soft and is varied by the prolate nature of the 
proton structure. The present results suggest the smaller neutron deformation /3„ than the proton deformation f3p as 
Pn < Pp in the ground band of ^"'Ne. On the other hand, /3„ is as large as (3p in the excited band. 

In order to extract information about the proton and neutron deformations from the experimental data, it is useful 
to analyze E2 transition strengths and compare them with those of mirror nuclei. Similar analyses on the electric 
moments were done in Ref. [20,21] where deformations of C isotopes have been discussed. Recently, B{E2; 2^ 0'^) 
in ^^Si has been measured by Coulomb excitation [4] as B(E2) = 4.6± 1.4 w.u., which is smaller than or of same order 
of B{E2; 2^ 0^^) =6.8± 1.6 w.u. in ^^Ne. If we assume the mirror symmetry between 24Ne and "^^Si, B{E2;'^'^S\) 
in ^^Si should equal to S„(S2;^^Ne) for neutron transitions in ^^Ne. Therefore, these facts suggest that, in ^^Ne, 
the neutron transition strength i?„(£'2;^^Nc) is as small as the proton transition strength _B(i?2;^"'Nc). The ratio 
B„(£;2;24Ne)/S(£'2;24Ne)< 1 is smaller than the naive expectation, {N/Zf « 2, given by a collective model picture. 
The reduction of B„(£;2;2^Ne) can be described by the smaller neutron deformation /3„ than the proton deformation 
fip in the 'oblate' ground band. In fact, the calculated ratio Bn{E2)/B{E2) = B{E2;'^^Si) / B{E2;'^'^Nc) in the 'oblate' 
band of ^^Ne is approximately 1 (table II), which well agrees to the experimental data. On the other hand, in 
the excited prolate band with /3„ w /3p, the calculated B{E2f^Si*) is twice larger than the B{E2f'^'Nc*), and the 
ratio Bn{E2)/B{E2) = B{E2-^'^^\*) / B{E2-^'^^e*) is consistent with the collective model expectation {N/Zf w 2. 
Comparing the calculations and the experimental data of B{E2-^'^'Hc) and _B(£^2;^^Si), we conclude that the 'oblate' 
band correspond to the ground band of ^''Ne. The present result of the 'oblate' band still overestimates the B{E2\'^'^S\) 
as shown in table II. It is conjectured that the neutron deformation might be smaller due to the d5/2-shell closure 
effect. 



D. Lowering mechanism of oblate states 

The shape coexistence in ^^Si has been studied by 7a-cluster models for a long time [7,8]. The oblate and prolate 
solutions coexist as local minima in the a-cluster model space, however, most of the calculations fail to reproduce the 
order of the oblate and prolate states in ^*Si. Even in the three-dimensional calculations of 7a clusters, the prolate 
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solution is lower than the oblate solution [8], which is inconsistent with the fact of the oblate ground state. Exceptions 
are the la calculations with a few parameter sets of Brink-Beoker forces in Rcf. [7], though the binding energy is 
much underestimated. 

The main reason for the failures in describing the deformation of the ground state of ^^Si in the 7a-cluster calcula- 
tions is because the spin-orbit force is omitted in the a-cluster models. As seen in the Nilsson orbits [3], the spin-orbit 
force plays an important role to make the shell gap at = 14 in the oblate system. It means that the spin-orbit force 
is significant to gain the energies of the oblate states in = 14 systems. Here we explain the lowering mechanism of 
the oblate state in ^^Si within the AMD model space. 

As described in VA, in the results of ^^Si it is found that the positions Di(thc real part of Xi/y^) of Gaussian 
centers correspond to the 7a-like cluster configurations (Fig. 6). If the imaginary parts of the Gaussian centers are 
zero, the expectation value of the spin-orbit force vanishes because it is exactly zero in the systems composed of the 
ideal (Os)'*-a clusters. However, in the intrinsic state of the oblate band, wc find that the non-zero imaginary parts 
of the Gaussian centers cause the rotational motion. The positions and the momenta of the singlc-nucleon Gaussian 
wave packets arc given by the real parts and imaginary parts \Kx\i of the Gaussian centers (Xi/^/z/) as given 

in Eq.5. {Dy\i and {Kx\i in the intrinsic state of the oblate ground band of ^*Si are plotted in Fig. 7. The positions 
and momenta for up-spin and down-spin nucleons are shown by square and circle points, respectively. The correlation 
between Dy and indicates the rotational motion around the ^;-axis. Namely, the angular momenta of the up-spin 
nucleons are parallel to the ^-axis and those of the down-spin nucleons are anti-parallel. Therefore, the rotational 
motions for up-spin nucleons and the down-spin nucleons are reverse to each other, and make the spin-orbit force 
attractive. The energy gain of the spin-orbit force is more than 20 McV as shown in table I. In the fc-th a-like 
cluster, the momenta K^iX) for the up-spin proton and neutron are opposite to those K^iX) for the down-spin proton 
and neutron as K^iX) = —J^xiX) ^ shown in Fig. 7. Although the positions Dj are located in the 7a-like cluster 
configurations, it is important that the a-likc clusters arc dissociated due to the imaginary parts K^, because the 
Gaussians for the up-spin nucleons are moving in the reverse to the down-spin nucleons in each of a-like clusters. 

Thus, in the AMD framework, the spin-orbit force can be gain while keeping the a-like cluster configurations. The 
energy gain of the spin-orbit force arises from the flexibility of the AMD wave functions where Gaussian centers are 
expressed by 'complex' variational parameters instead of real parameters. 

VI. SUMMARY 

We studied the deformations of iV = 14 isotones, ^^B, ^°C, ^^O, ^"'Ne and ^^Si, while focusing on the shape 
coexistence of the oblate and prolate neutron structure. The relative energies between the oblate and prolate states 
are found to be sensitive to the strength of spin-orbit force. By using a set of interaction parameters m = 0.62 and 
uis = 2800 MeV in the MVl(case 1)4-G3RS force, we can describe the properties of the shape coexistence in ^®Si, 
which is known to have the oblate ground band and the prolate excited band. The present results reasonably agree 
to the experimental data of radii, moments, and E2 transition strengths. 

In the present calculations, the solutions with the oblate neutron structure are energetically favored and form the 
ground bands in these A'' = 14 isotones. In the results of ^^B, the prolately deformed excited band with cluster 
structure was predicted to appear above the oblate ground state. In ^"^Nc, we proposed a possible difference between 
proton and neutron deformations in the ground band. We discussed Q-moments and El transition strengths in 
relation to the intrinsic deformations. In the analysis of the results, it was found that the spin-orbit force plays an 
important role in the lowering mechanism of the oblate neutron structure in A" = 14 systems. We found that the 
imaginary parts of the single-nucleon Gaussian centers in the AMD model space are important to gain the spin-orbit 
force. 
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FIG. 7. (a) Positions Dy and momenta Ky of 28 Gaussian centers in the intrinsic wave function of the oblate solution of 
^*Si. The X, y and z axis are chosen to be (x^) < {y^) < {z^)- It is found that the z-sx\s in the oblate solution almost equals 
to the direction of the intrinsic-spins of single-nucleon wave functions, (b) Schematic figure for the rotational motion given by 
the positions and momenta of Gaussian centers. 
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